Influence of rare regions on magnetic quantum phase transitions 
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The effects of quenched disorder on the critical properties of itinerant quantum magnets are con- 
sidered. Particular attention is paid to locally ordered rare regions that are formed in the presence 
of quenched disorder even when the bulk system is still in the nonmagnetic phase. It is shown that 
these local moments or instantons destroy the previously found critical fixed point in the case of 
antiferromagnets. In the case of itinerant ferromagnets, the critical behavior is unaffected by the 
rare regions due to an effective long-range interaction between the order parameter fluctuations. 
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Rare regions and their influence on observables is an 
important, if intricate, aspect of systems with quenched 
disorder. An effect that has been known for a long time 
is the formation of a Griffiths region [Q. To explain this, 
let us consider a ferromagnet for definiteness. Disorder 
will decrease the critical temperature from its clean value, 
T®, to a value T c < Xj? in the disordered system. In the 
temperature region T c < T < Tj? the system does not dis- 
play global order, but one will find regions that are devoid 
of any impurities and hence show local magnetic order. 
The probability of finding such a 'rare region' in general 
decreases exponentially with its size. The resulting mag- 
netization fluctuations have very slow dynamics. They 
are often called 'local moments' or 'instantons', and they 
lead to a nonanalytic free energy for all temperatures be- 
low T~, even though no long-range order develops until 
the temperature reaches T c . For generic classical sys- 
tems this is a weak effect, the singularity being only an 
essential one. An important exception is the model stud- 
ied by McCoy and Wu 0, which is a two-dimensional 
(2-d) Ising model with random bonds in one direction, 
but identical bonds in the second direction. The infinite 
correlation of the disorder in this model leads to much 
stronger effects, with the average magnetic susceptibility 
diverging in a finite- width temperature region above T c . 
The transition at T c is nevertheless sharp. The diver- 
gence of the average susceptibility for T > T c is caused 
by atypical fluctuations in the susceptibility distribution, 
and the averaged order parameter becomes nonzero only 
for T < T c . The temperature region T c < T < T c ° is 
known as a Griffiths region. Little is known about the 
influence of rare regions on the critical behavior at T c , 
and in the conventional theory of the critical behavior of 
disordered magnets jj| the rare regions are neglected. 

Recent work M on a random- T c classical Ising model 
has suggested that the effects of the rare regions go be- 
yond the formation of a Griffiths region, even in this sim- 
ple model where the conventional theory || predicts stan- 
dard power-law critical behavior. These authors showed 
that the conventional theory is unstable with respect to 



perturbations that break the replica symmetry. By ap- 
proximately taking into account the rare regions, they 
found a new term in the action that actually induces 
such perturbations. In some systems replica symmetry 
breaking is believed to be associated with activated, i.e. 
non-power law, critical behavior. Although no final con- 
clusion about the fate of the transition could be reached, 
Ref. thus raised the possibility that the random- T c 
classical Ising model shows activated critical behavior as 
a result of rare-region effects, as is believed to be the case 
for the random-field classical Ising model |J. 

The problem of rare regions is even less well inves- 
tigated for the case of quantum phase transitions, i.e. 
transitions that occur at T — as a function of some 
non-thermal control parameter . An important excep- 
tion to this are certain l-d systems. Fisher M has inves- 
tigated quantum Ising spin chains in a transverse random 
field, which is closely related to the classical 2-d McCoy- 
Wu model (with time playing the role of the dimension 
along which the disorder is correlated). He found ac- 
tivated critical behavior due to rare regions, which has 
been confirmed by numerical simulations § . Recent sim- 
ulations J^] suggest that this type of behavior may not 
be restricted to l-d systems, raising the possibility that 
exotic critical behavior dominated by rare regions may 
be generic in quenched disordered quantum systems. 

Apart from their relevance for disordered magnets and 
their critical properties, rare regions are believed to be a 
crucial ingredient for understanding other systems with 
quenched disorder. For instance, it has been proposed 
that a complete understanding of the properties of doped 
semiconductors, and of the metal-insulator transitions 
that are observed in such systems as a function of doping, 
requires the consideration of local moments |[o]-[l2| . 

In this Letter we study this important problem analyt- 
ically for quantum phase transitions in d > 1. We con- 
centrate on magnetic transitions, and contrast the cases 
of itinerant ferromagnets (FMs) and antiferromagnets 
(AFMs), respectively. We find that for the latter, the 
rare regions destroy the critical fixed point (FP) found 
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in a previous study , and thus have a profound effect 
on the critical behavior. In contrast, for itinerant FMs we 
find that, for certain realizations of the disorder, the pre- 
viously found critical behavior [|l4| is stable with respect 
to rare regions, due to an effective long-range interaction 
between the spin fluctuations. In addition, we find that 
the ultimate effects of the rare regions depend on how the 
disorder is realized in a particular system. Therefore, no 
generally valid conclusions are possible, and the effects 
of rare regions must be studied carefully and separately 
for each system under consideration. 

Let us first consider the case of itinerant quantum an- 
tiferromagnets. Our starting point is the same as in 
Ref. fill , namely Hertz's action, which is a 4 -theory 
for a p-component order parameter field cf> whose expec- 
tation value is proportional to the staggered magnetiza- 
tion. The action reads 

S= dx dy <p(x)T(x,y) 4>{y) + u I dx (<p(x) ■ <p(x)) 2 . 



(la) 

Here x = (x, t) comprises position x and imaginary time 
r, and J dx = J dx dr. We use units such that 
h = Jcb = 1. r(x, y) is the bare two-point vertex function, 
whose Fourier transform is 



r(q,w n ) = (t + q 2 + K|)/2 



(lb) 



Here t denotes the distance from the critical point, q 
is the wavevector, u„ is a bosonic Matsubara frequency, 
and we measure both q and w n in suitable units. 

Disorder is introduced by making t a random function 
of position, t = to + St (it), where St(x) obeys a Gaus- 
sian distribution with zero mean and variance A. The 
standard procedure is to integrate out the 'random mass' 
5t(x) by means of the replica trick ||, which produces a 
term of order </> 4 with coupling constant A, in addition to 
the ordinary quantum fluctuation term in Eq. (^) with 
coupling constant u. The resulting theory does not easily 
allow for saddle-point solutions that are inhomogeneous 
in space, and to incorporate rare regions into it would be 
very difficult. We therefore follow a different procedure. 
In analogy to Ref. ^ , we consider inhomogeneous saddle- 
point solutions of the theory for a fixed realization of the 
disorder. The inhomogeneity comes about since St(~x.) 
has 'troughs' that make t < in some region in space, 
even though to > 0. Troughs that are sufficiently deep 
and wide support locally nonzero saddle-point solutions. 
These regions we will refer to as 'islands'. Outside of the 
islands, the solution is exponentially small. This means 
that for a system with N islands, and in the case of an 
Ising model (p = 1), there will be 2 N almost degenerate 
saddle-point solutions that can be constructed by con- 
sidering all possible distributions of the sign of the order 
parameter on the islands. For p > 1 there is a whole 
manifold of almost degenerate saddle points. 



Let <&(x) be one of these saddle-point solutions, and 
let us consider fluctuations about it, 4>(x) = «&(x) + tp(x) 
p5j . The different saddle points are far apart in config- 
uration space and separated by large energy barriers. If 
we restrict ourselves to small fluctuations about each sad- 
dle point, we can therefore write the partition function 
approximately as the sum of all contributions obtained 
from the vicinity of each saddle point, 

j D[*(x)] P[*(x)] J D[<p(x)]e- s l*W + VWl, (2) 

where P[<fr(x)] denotes the distribution of saddle points, 
and J indicates an integration over small fluctuations 
only. It is clear that this approximation takes into ac- 
count effects that one would call 'non-perturbative' in a 
standard treatment of quenched disorder. Also, consis- 
tent with our approximations, it can be shown that the 
inhomogeneous saddle-point solutions lead to a lower free 
energy than the homogeneous saddle point <j>(x) = 0. 

Performing the integration over the in Eq. (g) explic- 
itly is very difficult, and the result will in general depend 
on the properties of the distribution function P, which 
in turn depend on the details of the miscroscopic dis- 
order. However, a very basic observation simplifies our 
task: The 3?(x) represent static randomness, and the av- 
erage over this randomness is performed for the partition 
function. That is, we are dealing with static, annealed 
disorder. This is physically sensible, as the local moments 
are a self-generated part of the system and therefore in 
equilibrium with the rest of the degrees of freedom [ p"5| . 
In addition, of course, there is quenched disorder due to 
the underlying random mass term. This we handle by 
means of the replica trick. If we assume that the distri- 
bution of the 4> is short-range correlated (which will be 
the case for certain classes of realizations of the disorder, 
but not for others), we can immediately write down the 
effective action up to and including terms of 0(ip 4 ): 

S cS = ]T f dxdy ip<*(x)T Q (x,y)<p a (y) 
+u^2 I dx (<p a (x) ■ V a (x)f 

a •* 

-^(A + ™<5 a/3 ) fdxdy5(x-y) ( V a (x)) 2 

x(^(y)) 2 + 0( V 6 ) . (3) 



Here Tq is the Gaussian vertex, Eq. (lb), with t = to, 
A is the variance of the Gaussian random mass distribu- 
tion, and a and (3 are replica indices, w is the coupling 
constant of the annealed disorder term. We have also 
derived Eq. (||) by means of a detailed technical proce- 
dure which will be reported elsewhere Eq] . The technical 



derivation shows that w has the form w 



' v, with v 



characteristic of the distribution P, and it also yields 
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terms of 0((f 6 ) and higher. These turn out to be less 
relevant for the critical behavior than the quartic terms 
shown in Eq. (|J) . Notice that the annealed disorder con- 
tribution becomes indistinguishable from the usual cp 4 or 
it-term in the case of a classical transition. This is the 
reason why the authors of Ref. Q , who studied classical 
magnets, considered replica symmetry breaking in order 
to describe nontrivial effects of the rare regions. In the 
quantum case we get a nontrivial effect even at the level 
of a replica symmetric theory, which means that the in- 
fluence of rare regions on quantum transitions is stronger. 

To discuss the properties of the effective action, Eq. 
(||), we proceed as in Ref. jljj. We consider d = 4 — e 
space dimensions and e T time dimensions, and control 
perturbation theory by means of a double expansion in 
e and e T fT?j| . Defining w — wT~ eT , and putting T = 0, 
we obtain the following renormalization group (RG) flow 
equations to one-loop order, 

^ = (e - 2e T )u - 4(p + 8)u 2 + 48uA , (4a) 
dl 

^ = eA + 32A 2 -8(p + 2)uA + 8pAw , (4b) 

^ = (e - 2e T )w + Apw 2 - 8{p + 2)uw + A8Aw . (4c) 

An analysis of Eqs. (^) shows that they possess eight FPs. 
Four of them have a vanishing FP value of w, w* = 0, 
and have been discussed before in Ref. . Of particular 
interest is the nontrivial critical FP u* — (e + e T )/16(p — 
1), A* = [(4-p)e + 4(p + 2)e T ]/64(p-l), w* = 0, which 
on the w = hypersurface is stable for p smaller than 
some p c . To one-loop order, and for e = e T , p c = 16. A 
linear stability analysis reveals that the third eigenvalue, 
\ a = (4 - p)(e + 4e T )/4(p - 1), is positive for p < 4. 
In the most interesting case p = 3, w is thus a relevant 
operator with respect to this FP, which means that the 
rare regions destroy the FP. It is, however, interesting 
to note that for p > 4 the FP is stable and describes 
power-law critical behavior. There also are four FPs with 
w* =/= 0. They are all unstable except for one with w* = 
(p — 4)(e + 4e r )/8p(10 — p), which is negative for p < 4. 
Since the bare value of w is positive, and the structure of 
the flow equations does not allow for w to change sign, 
this FP is unphysical. There is thus no new FP for p < 4, 
and a numerical solution of the flow equations reveals 
runaway flow in all of physical parameter space. 

We conclude that for p < 4 the AFM long-range or- 
der found in Ref. |L3) is unstable against effects induced 
by rare regions, a result that is consistent with the pre- 
vious suggestion that AFM long-range order is strongly 
suppressed by quenched disorder Jl0| . However, other 
possibilities exist. For instance, there could be a tran- 
sition to a long-range ordered state, but with activated 
critical behavior which manifests itself as runaway flow 
in a perturbative RG calculation. The viability of this 



latter suggestion is underscored by the fact that a cal- 
culation of the local moment contribution to the order 
parameter susceptibility yields xlm(T) ~ 1/T |Q. This 
is similar to Fisher's 1-d result x(^) ~ l/^ 7 wfth 7 < 1 
M . (Our exponent value of unity is a result of our saddle- 
point approximation for the local moments.) This shows 
that we are really describing a Griffiths region, which was 
shown in Ref. to lead to a transition with activated 
critical behavior in d = 1. A third possibility is that a 
conventional critical FP exists, but cannot be described 
with perturbative RG methods. This possibility is con- 
sistent with the stability of conventional critical behavior 
against w for p > 4, as discussed above. 

We now turn to the case of itinerant ferromagnets, 
which constitute an interesting contrast to the AFM case. 
In Ref. jl4| it was shown that a description of itinerant 
FMs that neglects rare regions leads to an action that 
has the same form as Eq. (||) with w — 0, except that the 
bare two-point vertex function reads 



r FM 
1 



(q,w„) 



ld-2 



l/q 2 ) /2 



(5) 



and that the field ip(x) now describes ferromagnetic fluc- 
tuations. There are two crucial, and related, differences 
between Eq. (^) and its AFM counterpart. The first one 
is the structure of the frequency dependence, which en- 
ters as \u n \/q 2 [Q and reflects the diffusive nature of 
the spins in a disordered environment. In Ref. it 
was shown that the same diffusive spin dynamics leads 
to the |q| d ~ 2 term, which dominates the usual gradient 
squared term as long as d < 4. In the original treatment 
of quantum FMs by Hertz ]l8| ], loop corrections would 
have been required to find this term, while the method 
of Ref. |Q builds it into the bare theory. Consequently, 
the correlations between the spin density fluctuations are 
effectively long-ranged, a feature that is well known to 
stabilize the Gaussian critical behavior |L!|. Indeed, it 
was shown in Ref. |Q that the Gaussian critical behav- 
ior, with 77 = 4 — d, v = \/(d — 2), 7=1, and z — d, is 
stable for 2 < d < 4. Here 77, i/, and 7 are the usual crit- 
ical exponents, and z is the dynamical critical exponent. 
They can all be simply read off Eq. (Hl The exponents 
(3 and S were also determined in Ref. |14J] , their values in 
d = 3 are (3 — 2,6 — 3/2. The remarkable claim of Ref. 
Jlij was that these exponent values, which in d = 3 are 
very different from both mean-field values and classical 
Heisenberg values, constitute the exact critical behavior 
of itinerant quantum FMs. 

An obvious question is whether this claim survives the 
consideration of rare regions. To answer this, we perform 
an analysis analogous to the one for AFMs above. A 
simple way to incorporate the rare regions into the action 
is to write the quenched disorder or A-term in the action 
as a random mass in the Gaussian vertex (i.e., to 'undo' 
the integrating-out of the random mass), to construct 
inhomogeneous saddle-point solutions and expand about 
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them, and then to integrate over the manifold of saddle 
points as in the AFM case. Clearly, this leads to a w-teim 
in the action, like in Eq. (||). We have derived the same 
result starting from a more microscopic formulation. We 
will report the details of the derivation elsewhere p6[, 
here we mention only one important point: After Eq. (0) 
we mentioned that w is proportional to u 2 . Since u in the 
FM case is wavenumber dependent and diverges in the 
short- wavelength limit (i.e., its bare scale dimension is 
negative) this raises the question whether the bare 
value of w is finite. The answer is affirmative, since the 
ui-term arises from field configurations that are nonzero 
only on islands. The u that contributes to the bare value 
of w therefore has to be taken at wavenumbers that are 
on the order of a typical inverse island size, and hence 
is finite. Once again it is important here that the island 
size distribution falls off exponentially for large sizes. We 
can thus treat w as a number. 

Now let us perform a power counting analysis to de- 
termine the stability properties of the Gaussian FP. As- 
signing a length L a scale dimension [L] = — 1, the scale 
dimension of the imaginary time is [r] = — z — — d. For 
the scale dimension of the field we find [<p(x)] = (d+2)/2. 
The scale dimensions of both w and A then become 
[w] = [A] = 4 — d, i.e. they are irrelevant with respect 
to the Gaussian FP for d < 4, and marginal in d = 4. 
Terms of higher than quartic order in cp that are pro- 
duced by a technical derivation of the effective action 
p6[ turn out to also be irrelevant. We thus conclude 
that the FM critical behavior determined previously |l4| ] 
is stable against rare regions physics, in sharp contrast to 
the AFM case. The reason for this qualitative difference 
is the effective long-ranged interaction between the order 
parameter fluctuations (as expressed in Eq. (^) and in 
the value of the exponent 77), which is sufficient to sup- 
press all disorder fluctuations, including the ones due to 
rare regions. By the same arguments, the FM Gaussian 
FP is also stable against replica symmetry breaking. 

We conclude with one additional remark. One might 
ask why the rare regions or local moments don't cut off 
the singular wavenumber dependences |q| d ~ 2 and |o; n |/q 2 
in the Gaussian vertex, Eq. (H) . The reason why this does 
not happen is that both singularities are consequences of 
spin diffusion, which in turn is a consequence of the spin 
conservation law. The rare regions ultimately derive from 
a spin-independent disorder potential, which clearly can- 
not destroy spin conservation. We note, however, that 
the above arguments are restricted to a tree-level analysis 
of our effective field theory. Although the effective theory 
is a sophisticated one, which at tree level contains many 
effects that would require loops in more standard treat- 
ments, we of course cannot exclude the possibility that 
loop corrections might lead to qualitatively new terms in 
the action. If such new terms included a RG-generated 
spin dependent potential, then this might change our con- 
clusions. However, at such a level of analysis one would 



also have to include effects due to interactions between 
the rare regions, which we have mostly neglected. Such 
interactions are known to weaken the effects of the rare 
region , but in general it is not known by how much. 
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